ADP-03-109/T547 



Delta Baryon Magnetic Moments From Lattice QCD 
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Theoretical predictions for the magnetic moments of the physical A baryons are extracted from 
lattice QCD calculations. We utilize finite-range regulated effective field theory that is constructed 
to have the correct Dirac moment mass dependence in the region where the u and d quark masses 
are heavy. Of particular interest is the chiral nonanalytic behaviour encountered as the vr decay 
channel opens. We find a A"^"^ magnetic moment (at the A pole) of = 4.99 ± 0.56 /xjv. 

This result is within the Particle Data Group range of 3.7-7.5 /xjv and compares well with the 
experimental result of Bosshard et al. of /iA++ = 4.52 ± 0.51 ± 0.45 hn- The interplay between the 
different pion-loop contributions to the A+ magnetic moment leads to the surprising result that the 
proton moment may exceed that of the A"*", contrary to conventional expectations. 

PACS numbers: 12.39.Fe, 12.38.Gc, 13.40.Em, 14.20.Gk 



I. INTRODUCTION 

The magnetic moments of the A baryons have already 
caught the attention of the experimental community and 
hold the promise of being accurately measured in the 
foreseeable future. Experimental estimates exist for the 
A++, based on the reaction 7r+ p — » 7r+ 7' p. The Par- 
ticle Data Group ^Ij] cites a range of values, 3.7-7.5 hm, 
for the A++ magnetic moment, with the two most recent 
experimental results being = 4.52 ± 0.51 ± 0.45 /ijv 

and ^A++ = 6.14 ± 0.51 This gives an idea of 

the elusive nature of the A magnetic moments. In princi- 
ple, the A+ magnetic moment can be obtained from the 
reaction 7 » ^ 7r° 7' p, as demonstrated at the Mainz mi- 
crotron [J, 13 and KotuUa et al. have recently reported an 
initial measurement of /iA+ = 2.7l|^'3(stat.)±1.5(syst.)± 
3(thcor.) 0- 

Recent extrapolations of octet baryon mametic mo- 
ments have utilized a Pade approximant 0,11, 0, Il0l |. 
an analytic continuation of chiral perturbation theory 
(xPT) which incorporates the correct leading nonana- 
lytic (LNA) structure of xPT. This idea was recently 
extended to a study of decuplet baryons in the Access 
Quark Model (3> 

where the next-to- leading nonanalytic 
(NLNA) structure of xPT was also included. Incorpo- 
rating the NLNA terms into the extrapolation function 
contributes little to the octet baryon magnetic moments, 
however it proves vital for the decuplet. This is because 
the NLNA terms from xPT contain information regard- 
ing the branch point at the octet-decuplet mass-splitting, 
TOtt = 5, associated with the A Nn decay channel, 
which plays a significant role in decuplet magnetic mo- 
ments. Here we include the NLNA behaviour in a chiral 
extrapolation of the magnetic moments of the A baryons 



calculated in lattice QCD. 

Although the Pade technique provides a well behaved 
method of extending to heavy quark masses, it is 

not so obvious how it could be extended to higher order 
in a chiral expansion. We have therefore chosen to cal- 
culate explicitly the pion loop diagrams which give rise 
to the LNA and the NLNA behaviour of the A rnagnetic 
moments using finite-range regularization p^ . Il3j . In 
this investigation we select the dipole-vertex regulator. 
We stress that, in the context of the chiral extrapolation 
of lattice QCD data, the use of a finite-range regulator 
is the preferred alternative to dimensional regularization 
as demonstrated in Ref. Ideally, the regulator mass 
should be constrained by lattice QCD results. However, 
the data available lies at large quark masses and pro- 
vides little guidance. As such, we explicitly explore the 
regulator-mass dependence of our results, finding a sensi- 
tivity of 3.6% or less in the charged A magnetic moments, 
for variation of the dipole mass parameter in the range 
0.5-1.0 GeV. 



II. EXTRAPOLATIONS - AN EFFECTIVE 
FIELD THEORY 



To one-loop order the dimensionally-regulated yPT re- 
sult for decuplet baryon magnetic moments is [Tj, [l5| 



M 



N 



+ T{-5, TOj, /^)) +a2ml. (1) 
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The J- {5, nij , ^) functions (with S the octet-decuplet mass 
difference) are the nonanalytic components of the meson 
loop diagrams shown in Fig. 1. These expressions are 
given explicitly in Rcf. 14] . The constants oq and 02 
multiplying the analytic terms are not specified within 
xPT and must be determined by other means. 
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TABLE I: The chiral coefBcients for each intermediate state 
for A baryon transitions. 
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FIG. 1: Diagrams indicating the leading and next-to-leading 
contributions included in the calculations. The letters T 
and B correspond to decuplet and octet baryons respectively. 
The dashed line is the appropriate meson accompanying the 
baryon transition. The analytic expressions for these dia- 
grams are given explicitly in Eqs. Q. 



It is well known that dimensionally-regulated xPT ex- 
pansions for baryon properties are troubled by a lack of 
clear convergence 0, 0, 0, This problem must 

be addressed if xPT results are to be apphed far beyond 
the chiral hmit, for example in the extrapolation of lattice 
data where pion masses up to 1 GeV are involved. In- 
deed, truncated dimensionally-regulated xPT expansions 
for magnetic moments exhibit the incorrect behaviour for 
large ttIt^. To address these concerns we adopt finite- 
range regulated xPT by calculating the meson-loop di- 
agrams using a dipole- vertex regulator. The loop con- 
tributions then approach zero naturally as becomes 
large. We combine the analytic terms of the chiral ex- 
pansion into a single term that maintains the chiral ex- 
pansion to the order we are working, while guaranteeing 
the correct magnetic moment mass dependence at heavy 
pion masses. 

The extrapolation function which wc employ to obtain 
decuplet baryon magnetic moment predictions from lat- 
tice QCD data, thus takes the form 



a-o 



I + Co ml 



^Xt' GTT'{mj,A) 

T' 

+ Y.^BGTB{m„6,A), (2) 



where the sums over T' and B include the decuplet and 
octet intermediate states, respectively. The first term en- 
capsulates the pion mass dependence of the photon field 



coupling directly to the baryon and reflects previous suc- 
cess in the extrapolation of octet baryon magnetic mo- 
ments 013, El- This term ensures the correct behaviour, 
as a function of quark mass, of the Dirac moments of the 
quarks at heavy quark masses (as oc niq up to 1 GeV^ 
- above this range one should use niq directly). The pa- 
rameters ao and cq are chosen to optimize the fit to lattice 
data. The functions GTT'iTnj, A) and GTsijnj , A) are 
the (heavy baryon) loop contributions to the magnetic 
form factors, in the limit q ^ 0, for decuplet-decuplet 
and decuplet-octet transitions respectively, see Fig. 1. 
These functions are given by 

GTT'(m„ A) = lim — / d^fc U{k) U{k') /'^ ^ , 
9^0 27r J [ujk Wfc' y 

GTB{mj,S,A) = lim -!- / d^fc Uik) U{k') 
q^o 2n J 

[uok +i^k' - 5){q X kf 



uik (jJk' {uik + Wfe' ) (wfc - (5) {ujk' - (5) ' 

(3) 

where q, k, k' are the momenta of the photon, incoming 
meson and outgoing meson, respectively. Note, k' ~ k + q 
and ujk = ^k'^ + . The functions U{k) and U{k') are 

used to regulate the loop integrals - we use a dipole with 
a finite mass parameter A. Therefore 



U{k) = 
U{k') = 



A2 



A^ + k\ 

A^ ' 
A2 + {k') 



(4) 
(5) 



The masses, rrij, of Eq. jSJ are octet-meson masses as- 
sociated with the meson cloud and Xi &re model indepen- 
dent constants |23| that give the magnitude of the loop 
contribution for each intermediate baryon state. These 
coefficients are given by 14j 
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where H describes meson couplings to decuplet baryons 
and C is the octet-decuplet couphng constant. We assign 
n and C their 311(6) values of H = -3D and C -2D 
PH where the tree level value for D is 0.76. The decay 
constants take the values /tt = 93 MeV and Jk — 1.2 /tt 
[2^ . appropriate to an expansion about the chiral SU(2) 
limit. Note also that Mjv is the nucleon mass and the 
parameters a and /3 are given in Ref. p^ . The model 
independent loop coefRcients, x, are summarized in Ta- 
ble □ 

In these calculations the complex mass scheme 0, 01 
is adopted as a method of incorporating the finite life of 
the A resonances, whilst also retaining electromagnetic 
gauge invariance. That is, one aims to extract the mag- 
netic moment of the A at the pole position in the complex 
energy plane, (5^?°^°) = Sr - i 61 = 270 - i 50 MeV Q for 
the physical pion mass (mP^^*^). Since the value at the 
position of the A pole is independent of the path chosen, 
we illustrate extrapolations along the path constant 
and Sj given by 



A2 



A2 + 5| 



Q{6r - m^) . 

(7) 

The latter is motivated by the usual expression for the 
A ^ TT self-energy i with G chosen to ensure Sj = 
50 MeV at mPi^y". 

To relate the kaon and pion masses we utilize the fol- 
lowing relations provided by xPT 



m 
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4' 



1 



phys\2 



(8) 
(9) 



This allows us to also incorporate kaon loops (at fixed 
strange quark mass) in extrapolations as a function of 
m^, to the physical mass regime. 

If we expand the first term of Eq. |(21) as a Taylor se- 
ries about TOtt = to order O(m^), and take the limit as 
A — > 00 in the loop integrals, we obtain the traditional 
xPT expansion for decuplet magnetic moments, as given 
in Eq. (where the infinite constants encountered in 
dimensional regularization simply redefine oq). Further, 
since the loop contributions approach zero much faster 
than for any reasonable value of A, Eq. (|2Jl guar- 

antees the correct mass dependence of the Dirac magnetic 
moment in the heavy quark mass regime. From previous 
studies of the nucleon axial form factor 25] it has been 
consistently demonstrated that A in the dipole regulator 
must have a magnitude < 1 GeV. For this investigation 
we assign A = 0.8 GeV, the optimal value obtained in 
analyses of state of the art full QCD simulations of the 
nucleon mass [T^. 

In Fig. 2 we plot Eqs. (0) (without the chiral coefficient 
prefactors) with the above value of S and a regulator of 
A = 0.8 GeV. Four types of intermediate baryon states 
are considered. The opening of the N tt decay channel has 
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FIG. 2: Plots of the four types of loop contributions to the 
A magnetic moments (without the coefficients) that are con- 
sidered here. These plots are obtained from Eqs. where 
we have used the complex mass scheme as a gauge invariant 
method with which to incorporate the finite life-time of the 
A resonances. 
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FIG. 3: Plots of the total loop contribution considered for 
each A baryon. These curves are the sum of the functions 
plotted in Fig. 2 multiplied by the appropriate coefficients 
listed in Table |l| 



an interesting effect on the magnetic moment contribu- 
tion. Fig. 3 presents a plot of the total loop contribution 
for each A baryon, determined using Eqs. lO and the 
chiral coefficients given in Tabled 



III. RESULTS 

The lattice QCD results given in Table IH] are ex- 
tracted from Ref. (2^] . These lattice calculations employ 
sequential-source three-point function based techniques 
|27j utilizing the conserved vector current such that no 
renormalization is required in relating the lattice results 
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kappa 


Baryon Mass 


Pion Mass 


A++ 


A+ 


A- 


A- 




(GeV) 


(GeV) 


(Miv) 


(mat) 




(Miv) 


0.152 


1.74 (50) 


0.964 (12) 


2.81 (18) 


1.40 (9) 


0.000 (00) 


-1.40 (9) 


0.154 


1.57 (60) 


0.820 (11) 


3.19 (28) 


1.59 (14) 


0.000 (00) 


-1.59 (14) 


0.156 


1.39 (80) 


0.665 (13) 


3.67 (43) 


1.83 (21) 


0.000 (00) 


-1.83 (21) 



TABLE II: Quenched lattice QCD magnetic moments for the A baryons extracted from Ref. (2^. Uncertainties are obtained 
using a third-order single-elimination jack-knife error analysis. Lattice results for the A" are exactly zero, in the quenched 
approximation. 



to the continuum. These simulations utihzed twenty- 
eight quenched gauge configurations on a 24 x 12 x 12 x 24 
periodic lattice at /3 = 5.9, corresponding to a lattice 
spacing of 0.128(10) fm. Moments are obtained from 
the form factors at 0.16 GeV^ by assuming equivalent 
dependence for both the electric and magnetic form 
factors. This assumption will be tested in future stud- 
ies incorporating the ideas presented here. Uncertainties 
are statistical in origin and are estimated by a third-order 
single-elimination jack-knife analysis (28i] . 

In Figs. 4-7 we present fits of the extrapolation func- 
tion, Eq. ((SJ, to the lattice data as a function of mj. The 
resulting magnetic moment predictions, along with the fit 
parameters gq and cq, are summarized in Table ITTll In 
Fig. 4 we show two experimental values for the mag- 
netic moment, the result fJ.A++ = 4.52 ± 0.50 ± 0.45 
is from Ref. and )Ia++ — 6.14 ± 0.51 fiN is given in 
Ref. 0. The discrepancy between these two results is 
a reasonable indication of the current level of systematic 
error in the experimental determination of the A"'"^ mag- 
netic moment. Our prediction of fJ, a+-i- — 4.99 ± 0.56 /xjv 
agrees well with the first experimental result, however it 
lies slightly below the range of the second. We note that 
the approach used in Ref. (Jj does not respect unitar- 
ity and when this shortcoming is addressed the authors 
find IJ.A++ = 6.01 ± 0.61 /LtAT 0], resulting in a some- 
what better agreement between our prediction and this 
experimental result. The result reported here for the A"*" 
magnetic moment of 2.49 ± 0.27 is in agreement with 
the initial measurement of KotuUa et al. namely 
Ha+ = 2.7l^;3(stat.) ± 1.5(syst.) ± 3(theor.)/ijv, however 
the experimental error at this time is still very large. 

To address the issue of regulator-mass dependence we 
vary A between 0.5 and 1.0 GeV, in each case readjusting 
oq and Co to fit the lattice data and find only a slight 
regulator-mass dependence of 3.6% or less, for each of the 
charged A baryons. Ultimately, lattice results at lighter 
quark masses will constrain A and therefore reduce this 
uncertainty. 

The inclusion of the A° is more for completeness rather 
than a firm result, as the lattice data in quenched QCD 
yields identically zero in this case. By varying the reg- 
ulator A between 0.5 and 1.0 GeV we find the moment 
remains positive, with the order of magnitude of our re- 
sult unchanged. We await new, unquenched lattice data 
in order to obtain a stronger prediction for the A" mag- 
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FIG. 4: The fit (solid line) to the lattice QCD data (solid 
circles) for the A"'"^ magnetic moment using the extrapola- 
tion function, Eq. (|5J . The dashed line is the plot of the first 
term of Eq. (|5J - and reflects the analytic terms of the chiral 
expansion. The solid star (*) indicates the prediction of the 
extrapolated lattice data. The two most recent experimen- 
tal results are indicated by the two asterisks (*), where the 
smaller prediction is taken from Ref. (3| and the larger from 
Ref. 1. 



netic moment. 

The interesting feature of the plots given in Figs. 4-7 
is the cusp at m,r = S which indicates the opening of the 
octet decay channel, A — > Ntt. The physics behind the 
cusp is intuitively revealed by the relation between the 
derivative with respect to of the magnetic moment 
and the derivative with respect to the momentum trans- 
fer, g^, provided by the pion propagator l/{q'^ + m'^) in 
the heavy baryon limit. Derivatives with respect to are 
proportional to the magnetic radius in the limit q^ —> 



' M. 



= -6 



dGM{q^) 



dq^ 



9^=0 



(10) 



If we consider for example A^'^ — > p'K'^ ^ with \j^mj) — 
|3/2,3/2), the N it state is in relative P-wave orbital an- 
gular momentum with \l,mi) = |1,1). Thus the pion 
makes a positive contribution to the magnetic moment. 
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Baryon ao co Lattice (pn) Experiment (/ijv) 

A++ 4.87 0.82 4.99 (56) 3.7-7.5 

A+ 2.44 0.83 2.49 (27) 2.7ti !^(stat.) ± 1.5(syst.) ± 3(theor.) 

A'' 0.63 381 0.06 (00) 

A" -2.40 0.80 -2.45 (27) 



TABLE IIL Theoretical predictions for the four A baryon magnetic moments from extrapolated lattice QCD results. The 
quoted errors are statistical in origin. The fit parameters oo and cq are given for each scenario. The experimental values for 
the A^"*" and A"*" magnetic moments are taken from Ref. Q. and Ref. respectively. 




FIG. 5: The fit (sohd fine) to the lattice QCD data (solid 
circles) for the A"*" magnetic moment using the extrapolation 
function, Eq. l(2j . The theoretical prediction for the A""" mag- 
netic moment is indicated by a solid star (★). The initial 
experimental value for the A""" magnetic moment is discussed 
in the text. The proton extrapolation (dotted line) of the 
lattice QCD simulation results of Ref. |29j| (open circles) is 
included for reference. 



FIG. 6: The fit to the lattice QCD data for the magnetic 
moment of the A", using the extrapolation function, Eq. ||5J. 
Symbols are as in Fig. 4. There is currently no experimental 
value for the A" magnetic moment. 
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As the opening of the p 7r+ decay channel is approached 
from the heavy quark-mass regime, the range of the pion 
cloud increases. Just above threshold the pion cloud 
extends towards infinity as AE 0, AEAt ^ h and 
the radius of the magnetic form factor diverges similarly, 
{d/dq^)GM oc {d/dm'^)GM —oo from above thresh- 
old. Below threshold, Gm becomes complex and the 
magnetic moment of the A is identified with the real part. 
The imaginary part describes the physics associated with 
photon-pion coupling in which the pion is subsequently 
observed as a decay product. Therefore decuplet to 
octet transitions enhance the A magnetic moments when 
the octet decay channel is closed. However, as the decay 
channel opens this physics serves to suppress the mag- 
netic moment as the chiral limit is approached. These 
transitions are energetically favourable, making them of 
paramount importance in determining the physical prop- 
erties of A baryons. 

The inclusion of octet-decuplet transitions in octet 
magnetic moment extrapolations is less important. Sig- 
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FIG. 7: The fit to the lattice QCD data for the magnetic 
moment of the A~, using the extrapolation function, Eq. (|5J. 
Symbols are as in Fig. 4. There is currently no experimental 
value for the A~ magnetic moment. 
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nificant energy must be borrowed from the vacuum for 
these transitions to occur. The formulation of effective 
field theory with a finite-range regulator naturally sup- 
presses transitions from ground state octet baryons to 
excited state baryons. Physically, the finite size of the 
meson-cloud source suppresses these high energy contri- 
butions. Hence any new curvature associated with octet- 
decuplet transitions for octet baryons is small 7] . Indeed, 
the inclusion of the p ^ A tt channels in the proton mo- 
ment extrapolation of Fig. 5 increases the predicted mo- 
ment by only 0.05 from 2.61 to 2.66 ^a?. For this 
reason octet to decuplet transitions have been omitted in 
other chiral extrapolation studies of the octet magnetic 
moments |§.lloj. 

In the simplest SU(6) quark model with m„ = 
the A"*" and proton moments are equal. However most 
quark models include spin dependent q — q interactions 
that enhance the A+ magnetic moment relative to the 
proton. This phenomenology is supported by the lattice 
results j26| at large quark masses. Consequeiitly, previ- 
ous linearly extrapolated lattice predictions [2g, using 
the same lattice data, suggested that the A"*" moment 
should be greater than that of the proton. We find evi- 
dence that supports a different conclusion, with predicted 
values for the A"*" and proton moments lying close at 
2.49(27) i-iN and 2.66(17) hn, respectively. The primary 
reason for this result is the interplay between the three 
different pion loop contributions; A+ A++7r^, n7r+ 
and A° 7r+. The transition to a A++ tt^ state largely can- 
cels the n n'^ contribution, which dominates the proton 
magnetic moment. This means that the A° 7r+ transi- 
tion is the main loop contribution to the A"*" magnetic 
moment, where the A"*" — > A" 7r+ coupling is weak, rela- 
tive to that of the p —>■ n tt"*" transition. 

The proton magnetic moment extrapolation is included 
in Fig. 5, and provides an illustration of the importance 
of incorporating known chiral physics in any extrapola- 
tion to the physical world. Thus, an experimental value 
for the A+ magnetic moment would offer very impor- 
tant insight into the role of spin dependent forces and 



chiral nonanalytic behaviour in the structure of baryon 
resonances. 



IV. CONCLUSION 

Finite-range regulated xPT has been applied to the ex- 
trapolation of lattice QCD results for the decuplet baryon 
magnetic moments. The magnetic moments of the four 
A baryons at the resonance pole are determined. Exper- 
imental values exist only for the A"*""*" magnetic moment. 
A result for the A+ should be forthcoming in the near 
future. The Particle Data Group gives a range of 3.7-7.5 
fiN for the A++ moment, with the two most recent ex- 
perimental results being /iA++ = 4.52 ± 0.51 ± 0.45 fj,N 
2] and /iA++ 6.14 ± 0.51 fiN 3]. Our lattice QCD 
prediction of /J.A++ = 4.99 ± 0.56 compares well with 
the first of these experimental results. An interesting re- 
sult from this investigation is the prediction that the A+ 
magnetic moment lies close to and probably below that 
of the proton moment. Therefore a new high precision 
experimental measurement of the A"*" moment would of- 
fer valuable insight into spin dependent forces and chiral 
nonanalytic behaviour of excited states. The lattice data 
used in these calculations is now over 10 years old, de- 
termined with pion masses greater than 600 MeV. The 
arrival of new lattice data at lower pion masses is there- 
fore eagerly anticipated, and will help constrain the fit 
parameters and associated statistical uncertainties. All 
of this should be forthcoming in the next few years and 
offers an excellent opportunity to test the predictions of 
QCD. 
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